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ABSTRACT 



Context. The density of the ambient medium where the supernova remnant evolves is a relevant parameter for its hydrodynamical 
evolution, for the mechanism of particle acceleration, and for the emission at TeV energies. 

Aims. Using XMM-Newton X-ray observations, we present a study of the ambient medium density of the historical supernova remnant 
SN 1006. 

Methods. We modelled the post-shock thermal emission to constrain the ambient medium density. Our study is focused on the North- 
West and the South-East rims of the remnant, where the thermal emission dominates. We used a plane-parallel shock plasma model 
plus another component for the ejecta that are not negligible in the regions of our study. The importance of the synchrotron component 
is also studied. In order to improve statistics, we combined several observations of the remnant. 

Results. The density found in the South-East rim is low, roughly 0.05 cm"^, and seems to be representative of the rest of the remnant. 
However, in the North-West rim (close to the bright optical filament), the density is significantly higher (about 0.15-0.25 cm"^^). This 
confirms a picture of SN 1006 evolving in a tenuous ambient medium, except in the North-West where the remnant has recently 
encountered a denser region. 

Conclusions. A density this low is compatible with the non-detection of the remnant by the HESS gamma-ray observatory. The lower 
density in the South-East implies a higher shock speed of 4900 km/s, higher than that of 2890 km/s measured in the North-West. This 
new estimate of the velocity could increase the maximum energy that accelerated particles can reach to energies of about 1 PeV. 
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1. Introduction 

The thousand-year-old Supernova Remnant (SN R) SN 1006 is 
a typ e la SNR located at a distance of 2.2 kpc (IWinkler et al.l 
l2003i) . The high Mach number coUisionless shock of the remnant 
heats and ionizes the interstellar medium (ISM). Furthermore, 
the blast wave is thought to be an efficient accelerator of cosmic- 
rays (CR). 

Evidence that electrons are accelerated at the shock is 
found at several wavelengths. The radio synchrotron emission 
(Reynolds & Gilmore 1986) arises from electrons accelerated to 
energies o f at least 1 GeV. The X-ray synchrotron emission in 
SN 1006 (lKovamaetal.lll995h comes from electrons of ~ 20 
TeV. Evidence concerning proton acceleration is more scarce. 
The best hope is to detect the gamma-ray emission of pion decay 
resulting from the interaction of accelerated protons (and heavier 
ions) with the thermal gas. However, the TeV emission can also 
result from Compton scattering of electrons at ^ 20 TeV on the 
CMB background. CuiTently there are only a few shell SNRs de- 
tected at those energies (Cas A, Ve la Junior, RX J1713.7-3946; 
lAharonianet'al. 2001. ,20051 l2006h and SN 1006 is not yet one 
of them. 

To interpret the TeV emission, it is important to estimate the 
gamma-ray flux from the tto mechanism. One of the most im- 
portant parameters for pr edicting this flux is the g as density in 
the y-ray emitting region. ^senofontov et al.l (l2005l) have shown 
that, in order to be consistent with the non-detection of SN 1006 
by the HESS experiment, the ambient medium density must be < 
0.1 cm"-'. The ambient medium density is also a key parameter 
for the hydrodynamical evolution of the remnant. As SN 1006 
results from a type la SN, the surrounding medium is expected 



to be unaltered by the progenitor. Moreover, the remnant is rel- 
atively isolated, as it is about 500 pc above the galactic plane. 
Thus it is a good case study for comparing observations with 
theoretical models of CR acceleration. 

In the literature, the gas density around SN 1006 ranges 
from 0.05 < mh ^ 1.0 c m~^. The different e stimates were 
based on radio HI study dPubrier et alj l2002h, UV spectr a 
jLaming et al.lll996t|Kq rreck et aLllIOOl [R aymond et al."200' 



and X-ray emission (Hamilton et al. 1986; Winkler & Lonj 



ll997l:lLonge"tan i2003). In SN 1006, while the limb-brightened 
synchrotron is dominant in the North-East and in the South- 
West, the thermal X-ray emission from the shocked gas domi- 
nates in the South-East (SE) and in the North-West (NW). It is 
possible to deduce the characteristics of this shocked gas, and in 
particular its density, with X-ray observ ations. In the NW , such 
a study has already been carried out by iLong et al] (l2003l) with 
Chandra. 

Using XMM-Newton, we present a study of the density in 
the South-East and in the North-West by modeling the thermal 
emission. 



2. Data processing 

The first observations of the remnant's So uth-East region by the 
XMM-Newton satellite jRothenflug et al.l 12004.) were not very 
deep. The effective exposure time in this region was 6 ks (ObslD 
0111090601) in comparison with 33 ks for the North-West 
(ObsID 0077340101). A second observation was carried out to 
provide better coverage of this region, on 2005 August 21 and 
22 with the medium filter for a duration of 35 ks. Unfortunately, 
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Table 1. XMM-Newton observations used in this paper 







MOS exp.(ks) 


Revolution.ObsId 


Observation Date 


Total 


Good 


0305.0111090601 (SE) 


2001 August 8 


15.77 


5.95 


0306.0077340101 (NW) 


2001 August 10 


64.85 


33.41 


1044.0306660101 (SE) 


2005 August 21 


35.15 


10.8 



the proton flare contamination was high and the remaining ex- 
posure time is only 1 1 ks. Note that for this new observation, 
the CCD number 6 (no longer working) and the CCD number 
4 (very noisy) of the camera MOS 1 were not used. All observa- 
tions used for our study are listed in Table[T]with their remaining 
exposure time after screening. We used only MOS data due to its 
better spectral resolution. All the spectra used are the average of 
MOSl and MOS2 detectors. 

The data were processed using the Science Analysis System 
(SAS version 6.5). All fitting was caiTied out using unbinned 
spectra with the C-statistics implemented in Xspec (vl2.2.1). 
Binning was used for graphical purposes only and fixed at 3cr 
for all spectra. We have fitted and plotted all the data from 0.3 
keVto lOkeV. 

To clean proton flare contamination from the event files, we 
used the MOS-FILTER process available in the Extended Source 
Analysis Software (XMM-ESAS). The foUowing method was 
used for background subtraction. To take into account the astro- 
physical background, we extract a spectrum from a region out- 
side of the SNR (regions A and B on Fig. [T]i; for statistical rea- 
sons, the area selected is as large as possible. The instrumental 
background spectrum i s derived from a 960 ks blank sky ob- 
servations compilation (ICarter & Readll2007h in the same detec- 
tor area as the source and renormalized in the 10-12 keV band 
(where astrophysical photons are negligible) over the full field 
of view. 

3. Spectral analysis 

To estimate the density of the ambient medium surrounding SN 
1006, we model the thermal X-ray emission at the shock. This 
gives us the post-shock density which is related to the gas density 
by the compression ratio, which is 4 if particle acceleration is 
negligible. 

Our study is focused on the SE and NW regions where ther- 
mal emission dominates. Figure [1] shows the extraction regions 
used for the spectra. In the NW, we used two region s: one is the 
region noted NW- 1 in the study of the remnant by iLong et akl 
( 1200 3^ with Chandra (to relate our study to what has been done 
previously) and the second one is the faint emission region and 
is noted NWf. Because the emission is faint in this region, we 
used a box as large as possible (1.2'x7.3'). The NW-1 region, 
in comparison, is much smaller (17"x3'). In the SE (region ESE 
and SSE), the width and length were fixed to 1' and 6.5' respec- 
tively, as a compromise between good statistics and being close 
to the shock. 

3.1. Method 

Our method is based on fitting the spectral data with thermal 
models and infering the post-shock density from the emission 
measure (EM) of the shocked ISM. Using Xspec, we get the 
parameter J rieriudl which is proportional to 

the EM. In this expression, Q. is the solid angle of emitting gas 
contained in our extraction region and dl the integration element 




aOO.a 30.0 04:00.0 



30.0 15:03:00.0 30.0 
Right ascension 



02:00.0 30.0 01:00.0 



Fig. 1. EPIC MOS plus PN image in the 0.5 to 2 keV. The scale 
is square root. The regions N W- 1 and NWf are focused respec- 
tively on the bright and faint emission. The boxes named A and 
B represent the regions used for the astrophysical background 
for the NW and the SE respectively. 



along the line of sight. We assume that the shocked ISM fills all 
the volume inside the SNR. This implies a filling factor of 1 . 

For a solar abundance plasma, the electron and hydrogen 
density are related by = 1 .2 1 xne. Therefore, for a given value 
of norm (obtained from the fit), we can solve the previous equa- 
tion if we know the density profile to compute the integral. We 
assumed a Sedov radial profile for that purpose. An important 
advantage of this method is that the hydrogen density is propor- 
tional to the square root of the ratio norm/ 1. Thus the density is 
weakly sensitive to variations of EM when using different ther- 
mal models, and to eiTors due to line of sight estimate. The eiTor 
due to a filling factor of 1 is also minimized for the same reasons. 

3.2. What kind of model can we use ? 

The pertinence of our estimate of the density will depend on 
the quality of the spectral model for the shocked ISM. As the 
ionization of the shocked ISM is out of equilibrium, we need 
to use non-equilibrium ionization (NEI) models. Moreover, the 
variation inside the SNR of the ionization timescale, defined as 
T - J n^dt, is important. Starting at zero at the shock, the ioniza- 
tion timescale increases toward the center of the remnant down 
to ~ 0.9 ^shock for a Sedov model and then decreases. Therefore 
we integrate cells of gas with different r along the line of sight. 

To model the emission, we used a simple model of plane- 
parallel shock (PSHOCK) with constant temperature, which al- 
lows the use of a distribution of t instead of a single ioniza- 
tion timescale like in a standard NEI model. A plane-parallel 
shock model is a good local approximation for radially small re- 
gions close to the shock, where t is still increasing inwards. The 
projection of the radial profiles onto the sky is not modelled by 
PSHOCK. 

To check the consistency of our model, we also determine 
the density by comparing the maximum ionization timescale t 
resulting from the fit with the ionization timescale computed us- 
ing Sedov solutions (refered to later as the t method). Sedov 
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Table 2. Best data fit parameters obtained for different models 
in region ESE. 









Shocked ISM 




Shocked ejecta 


Synchrotron 


Parameters 


C-stat 


kT (keV) 


T (cm^^^s) 


norm 


kT 


Cut-ofF frequency (Hz) 


One-component^ 


1266.1 


1.77 


3.8E9 


4.9E-4 






Two-component 


770.1 


1.75 


4.7E9 


2.25E-4 


0.52 




Three-component 


738.6 


0.85 


1.7E8 


0.4E-4 


0.62 


4.5E16 



" A simple NEI model was also tried but gave a worse fit. 



solutions have the advantage of taking into account the tempo- 
ral evolution of the density and the shock velocity. At a given 
relative angular distance fl/^shock, the only parameter required is 
the age of the remnant, which is well known. The method using 
the norm parameter, which will be refered to later as the EM 
method, is much more stable and will be used for the results. 
The T method is useful for checking the consistency between the 
models we use. In the following subsections, we use the region 
ESE as a test region to try different spectral models. 

3.3. Shocked ISM model 

As a first step, we assumed that synchrotron emission is negli- 
gible and that all of the emission comes from thermal emission 
of the shocked ISM. We therefore used a PSHOCK model with 
solar abundances (Lodders 2003]) a nd we fixed the int erstellar 
absorption to A^h = 7 x 10^° cm"^ jDubner et al.ll2002l) . In the 
NEI model (vl.l in Xspec) the oxygen He-K 6 and e lines are 
not included but they seem to be necessary to reproduce the large 
"shoulder" of the He-like oxygen line complex. To take them 
into account, we added one gaussian at 0.723 keV to represent 
the oxygen series He-K S and higher 

The spectral fit of region ESE is shown in Fig.|2](top panel) 
and the resulting best-fit parameters are given in Table |2] (see 
"one-component"). The fit is poor and we clearly see that some 
emission lines like magnesium (at 1.35 keV) and silicon (at 1.86 
keV) cannot be reproduced by this model. Such emission lines 
require oversolar abundances and must come from the ejecta. 

3.4. Shocked ISM with ejecta model 

To reproduce the spectrum correctly, we added another compo- 
nent coiTesponding to the ejecta emission (the same PSHOCK 
model as above but with non-solar abundances). 

To reduce the number of free parameters in this two- 
component model, we fixed most of the abundances of the ejecta. 
As there is strong evidence that SN 1006 is a Type la SN, we can 
obtain an estimate of the synthesized mass of elements and thus 
of the abundances in the ejecta through numerical simulations 
dlwamoto et al.l [19991) . Since there is neither hydrogen nor he- 
lium in type la explosions, the abundances are given relative to 
oxygen with respect to cosmic abundances from .Lodders ( 2003.) . 
We fixed the abundances of elements not constrained by emis- 
sion lines in the spectrum to the following values: 0=1.0 (by 
definition); C=1.0; N=0.1; Ai-2.0; Ca=3.0. Abundances from 
elements with prominent emission lines (Ne, Mg, Si) were let 
free. 

SN la are expected to produce an iron mass of 0.6 Mq, which 
is equivalent to a relative abundance of 20. This value is not com- 
patible with our data because the iron emission lines are very 
faint in the observed spectrum. The value preferred by the fit is 
compatible with zero. In practice, we fixed it to Fe=0.1. We can 



only see the emission from th e shocked ejecta, and the spatial 
distribution of iron is complex jWang & Chevalierll2001h . In our 
extraction region there appears to be little iron. 

The fit is much improved in comparison with the one- 
component model with AC-stat=495 (see Table |2] for C-stat and 
best-fit parameters). We note that the ejecta component is re- 
quired for a good fit at Mg and Si emission lines, but that it also 
contributes to the emission at low energy and in particular to the 
oxygen lines. In the first model, all the emission was attributed 
to the shocked ISM, whereas the emission of the ejecta is impor- 
tant, as seen in Fig|2] (middle panel). Using the EM method for 
the shocked ISM component, we obtained a density «h=0.053 
(0.048-0.057) cm""*. The value derived from the t method is con- 
sistent with the EM method : nH=0.059 (0.05-0.07) cm^^ This 
model reproduces the spectrum well and the two methods give a 
consistent estimate of the density. 

3.5. Impact of a potential synchrotron component 

SN 1006 is well known to be a non-thermal X-rays emitter Most 
of the emission in the bright SW and NE limbs is due to syn- 
chrotron emission. In the regions of our study, the thermal emis- 
sion dominates at low energy, but can we completely neglect the 
non-thermal emission? 

In order to answer to that question, we built a three- 
component model by add ing a SRCUT mo del. We fixed the ra- 
dio spectral index to 0.6 dPver et al.ll200ll) . The radio flux was 
extracted from the same regions as the X-rays using the radio 
map from Rothenflug et al. (2004). This radio image is a combi- 
nation of VLA data at the average frequency of 15 17.5 MHz and 
single-dish Parkes radiotelescope observations. 

The spectral fit is improved with a AC-stat=43. The result- 
ing cut-off frequency of 4.5 x 10^^ Hz (about 200 eV) is co m- 
patible with the frequency found bvl Rothenflug et all (|2004|) in 
their Fig. 7 (at a position angle 130 degrees). As we can see in 
Fig|2](bottom panel), the synchrotron component dominates the 
spectrum at high energy, at low energy, the shocked ejecta com- 
ponent dominates. 

In the two-component model, the shocked ISM component 
was constrained by the spectrum at high energy, while the ejecta 
were constrained mainly by the Mg and Si emission lines. Here, 
as the synchrotron reproduces the high energy, the shocked ISM 
is no longer well constrained. In order to correctly reproduce the 
data, we need a continuum and a thermal model for the emis- 
sion lines. The continuum can be handled either by a high tem- 
perature thermal model or by the synchrotron model. When we 
include both in the same fit, it is hard to disentangle them. The 
emission lines can be more easily reproduced by the shocked 
ejecta as several abundances are left unconstrained. Thus the 
shocked ISM component is almost not needed. The best-fit for 
a model with only synchrotron and shocked ejecta has a C-stat 
that is higher by only AC-stat=2. 
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Fig. 2. Best-fit spectrum from region ESE for the three models 
discussed. Top: shocked ISM. Middle: shocked ISM (dotted line) 
plus ejecta model (dashed line). Bottom: same model as above 
but with a synchrotron component (solid line). In the bottom 
panel, we note that the shape of the shocked ISM component 
is similar to the synchrotron component. The shocked ISM is al- 
most not needed, and the power law shape is due to a very low 
ionization timescale r ^ 10^ cm^-'s (no emission lines). For the 
middle and the bottom panel, the global fit (not shown) is a very 
good fit to the data. 



Including the synchrotron in the model could lower the esti- 
mate of the density. However, using hydrodynamical arguments 
we can set a lower limit on the density of the ISM, as explained 
later in Section l54l Therefore we think that the statistically pre- 
ferred model is not the best physical model. 

In conclusion, there are good physical reasons to include the 
synchrotron, but the data are not able to estimate its contribu- 



tion. We decided to use the simpler two-component model. The 
resulting densities should be viewed as upper limits. 

3.6. Using different sizes of extraction region 

The width of our extraction regions ranges from 17" for region 
NW-1 up to r for the regions ESE and SSE. We checked the 
impact of varying the size of the extraction region on the density 
estimate. In order to do so, we reduced the width of the region 
ESE to 17", keeping it immediately behind the shock. Using the 
two-component model described in Sect. 13.41 we derive densi- 
ties from the small ESE region with the EM and r methods of 
0.065 (0.055-0.075) cm"^^ and 0.06 (0.048-0.09) cm"^^ respec- 
tively. The density estimate is fully consistent with that deter- 
mined with a larger ESE region (see Sect 13.41 ). Therefore the 
density estimate is not sensitive to the size of the extraction re- 
gion. 

4. X-ray results 

4.1. General results 

Following the spectral analysis of region ESE (Sect. |3]l, we 
now use the same model for all the regions: shocked ISM plus 
shocked ejecta (as in Sect. 13.41 ). The results for the four regions 
(ESE, SSE, NWf, NW-1) studied here are listed in Table[3l 

For all the regions, the ejecta component has greatly im- 
proved the quality of the spectral fit in comparison with a one- 
component model. The emission lines are then well reproduced 
with our model as shown in Fig. [3] 

The density obtained with the EM method in the SE rim 
(ESE and SSE region) and in the NW faint region (NWf) is glob- 
ally low (hh^O. 049-0. 053 cm"^). Moreover the plasma parame- 
ters of the shocked ISM obtained from NWf are similar to those 
found in the SE. This is co nsistent wit h the smooth appearance 
of the Ha image ( Winkler et al.l 12003 ). We also note that the 
temperature of the shocked ISM is comparable in those three 
regions: kFisM^ 1-75-2.1 keV. This ISM temperature is much 
higher than the temperature of the shocked ejecta (0.5 keV). The 
shocked ISM component dominates the spectrum at high energy. 

As we can see on the spectra in Fig. [3] the ejecta component 
is important below 2 keV. Its contribution to the Mg, Si and Ne 
emission lines is obvious, and we can also see that part of the 
emission from the oxygen lines is produced by shocked ejecta. 
In the NW regions , an overabundanc e of Ne is required. This has 
also been noted bv lLong et al.l (l2003h in their spectral analysis of 
the bright filament. 

4.2. The bright filament 

The bright filament seen in X-rays in the NW is coincident with 
a bright Ha filament. A possible explanation is a "recent" in- 
teraction of the expanding SNR with a denser ISM. The sudden 
lower shock velocity flattens the shape of the blast wave in this 
region as it is seen in optical and X-rays. The thermal emission 
at the shock increases because the shock propagates a denser re- 
gion and thus creates a brighter filament than in the SE. The fact 
that the shock is entering in a denser re gion in the NW is par ticu- 
larly visible in Ha. The radio study bv lDubner et alj (l2002l) also 
detected a higher density in the NW. In Ha and in radio, we also 
see the faint structure of the region NWf. A possible interpreta- 
tion of this faint border is that the remnant is expanding in front 
of (or behind) the bright filament in a region of lower density 
similar to that of the SE. 
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Table 3. Best-fit parameters obtained with the shocked ISM plus shocked ejecta model (two-component) for the different regions. 
The gas density derived with the EM and the t methods is also given (see Sect. 13. It . Due to the faint emission in region NWf, the 
temperature and the abundances of the shocked ejecta component (parameters in square brackets) were fixed to the values of region 
NW-1. The lower and upper bounds are given in parentheses at a 90 % confidence level. * : The density estimate for region NW-1 
assumes constant density and velocity rather than Sedov profiles (see Sect. 14.21 ). 



Region 


NW-1 


NWf 


SSE 


ESE 


krejecta (keV) 


0.12(0.10-0.13) 


[0.12] 


0.55 (0.5-0.65) 


0.52 (0.49-0.57) 


Ne/0 


14 (8-22) 


[14] 


0.24 (0.1-0.3) 


0.22 (0.1-0.3) 


Mg/0 


49 (25-80) 


[49] 


6.0 (4.0-8.0) 


4.7 (4.0-5.5) 


Si/0 


21 (5-40) 


[21] 


19.0(17.0-23.0) 


17.0 (14.0-20.0) 


krisM (keV) 


1.45 (1.55-2.1) 


1.75 (1.5-2.1) 


2.1 (1.6-2.5) 


1.75 (141-2.15) 


T of sliocked ISM (10' cm-^ s) 


4.45 (3.5-5.6) 


24 (2.15-2.7) 


6.0 (5.4-6.5) 


4.7 (4.0-5.6) 


norm (10"** cm"^) 


0.55 (0.4-0.75) 


2.56 (2.2-2.9) 


1.88 (1.95-2.5) 


2.25 (1.8-2.6) 


C-statistics 










Two-component model 


454.8 


814.3 


750.1 


770.1 


One-component model 


646.6 


840.5 


1210.4 


1266.1 


Gas density (10~^ cm"'') 










T method 


15* (11.2-18.5) 


3.0 (2.6-3.3) 


7.5 (6.8-8.1) 


5.9 (5.0-7.0) 


EM method 




4.9 (4.5-5.2) 


5.2 (4.9-5.6) 


5.3 (4.8-5.7) 




0.5 1 2 5 0.5 1 2 

Energy (keV) Energy (keV) 



Fig. 3. Best spectral fit for the four regions with the double-component model (see Sect. 13.4b . The global fit is shown in black, 
the shocked ISM component in red (dominant at high energy), and the shocked ejecta in green. The 0.723 keV line (Sect. lTJt is 
included in the global fit. 



To derive the density from the plasma parameters, we used 
a simple Sedov model to compute t and EM and then compare 
them with the best-fit parameters (see Sect. 13. lb . However, in the 
NW-1 region it seems that part of the shock wave has recently 
encountered a denser region. For this reason, the geometry of 
the shell is too complex and probably not well represented by 



a Sedov density profile. The density derived with this method 
is probably not reliable in this particular region. Nevertheless, 
with a simpler model (assuming constant density and constant 
ve locit y) we can derive the density from the ionization timescale 
(ILong et al.l l2003l) . Indeed, we know th e shock speed in this 
region (2890 km/s, IWinkleret al.]|2003h and know that region 
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NW-1 has a radial extent of 17" (corresponding to an elapsed 
time of 240 years since the gas was shocked, at a distance of 
2.2 kpc). Thus for t=4.5 (3.5-5.5) xlO*^ cm"^ s, the density is 
«H - 0.15(0.11 - 0.18) cm"^. We note that krisM is lower in 
region NW-1 (1.45 keV) than in the SE and NWf (1.7-2.1 keV). 
This is consistent with the fact that the density in region NW-1 
is higher. 

In comparison , the study of the bright filament region by 
iLong et al] (l2003h with Chandra found a density of 0.25 cm"^. 
To see if there are differences between the data observed by 
XMM-Newton and Chandra, we tried the same si mple solar- 
abundance PSHOCK model (described in Sect. 3.2 of lLong et al.l 
^03) on our data. The absorption column density they use is 
higher than our value: 8.9 vs 7.0 xlO^" cm . We modelled the 
data with both values. The best-fit parameters are given in Table 
ID The model with low A^h is statistically preferred with a AC-stat 
of 35. 



Table 4. Comparison of the best-fit parameters obtained with 
XMM-Newton and Chandra data in the NW-1 region. 







kr^ (keV) 


r(10''cm-^' s) 






XMM Chandra 


XMM Chandra 




=8.9xl020cm-2 


0.6 0.6 


6.8 6.9 




=7.0xl020cm-2 


0.7 


6.5 



Both data give similar parameters and the fact that we are not 
using the same value for A^h does not have an important impact 
on the plasma parameters. Thus, the difference in the estimate 
of the density is not due to the data; it is mostly due to the dif- 
ference of models. When taking into account the ejecta using a 
two-component model, the density derived is slightly lower than 
when using a simple solar abundance PSHOCK model. 

The difference of density betwe en the NW-1 and othe r re- 
gions was alread y pointed out b y iKorreck et al] (|2004|) and 
iLong et all (l2003h . iKorreck et al.1 (|2004 estimated the density 
in the NE to be < 0.06 cm"^, which is compatible with what we 
found in the SE. 

We conclude that a low value of ~ 0.05 cm"-' is probably 
representative of the ambient medium around SN 1006, except 
in the bright NW filament where the density is higher (0. 15-0.25 
cm"^, depending on the model). 

5. Discussion 

5.1. Comparison with otiner observations 

By modelling the post-shock thermal emission of the remnant in 
the SE and in the NW, we found a density of ~ 0.05 cm"^. Many 
other studies (local and global) have been carried out to measure 
this density and provide estimates in the range of 0.05 cm"^ to 
1.0 cm^l 

Most of the local ones were focused on the bright Balmer 
line fil ament in the NW. Usi ng HST observations of this fil- 
ament, [Raymond et al.l ( |2007|) derived a density in the range 
of 0.25-0.4 cm~^ by modelling the Hot brightness behind the 
shock. By using the high spatial resolution X-ray observato- 
ries, some studies also est imated the density in the NW re- 
gion. IWinkler&Long|([T997 ). with ROSAT, used a model for the 
X-ray surface brightness profile behind the shock; they found 
a density of 1.0 cm However, using the Chandra X-ray 
Observatory, iLong et al.l (j2003) studied the same region and 
found a much lower value. They derived a density of 0.25 cm""* 



from the ionization timescale in a region right behind the shock. 
The spectral analysis of the same region using XMM-Newton ob- 
servations gave an identical value when using the same model 
(see Sect. 14. 2t . A model taking into account the ejecta compo- 
nent gave a slightly lower value (0.15 cm"^). We note that the 
density derived from the X-rays is lower than that derived from 
the optical. By using X-ray observations, we average the den- 
sity on the size of the extraction region, whereas in optical we 
measure the density right at the shock. If the shock wave is en- 
tering into a cloud, it is plausible that the difference of density 
between optical and X-rays is related to the gradient of density 
when entering the cloud. 

There is a general agreement that the remnant is encounter- 
ing a denser region in the NW, so the density there is not directly 
comparable with what we found in the SE. 

The global studies average the density all around the rem- 
nant. A surrounding density of ~ 0.3 cm~^ has been inferred by 
iDubner et al.l(l2002l) with a survey of the HI emission. However, 
this result is the average density along a line of sight of several 
hundreds of pc, which is very large in comparison with the typi- 
cal size of the remnant of 20 pc. 

[Hamilton et al.l (Il986l) derived a density of 0.05 cm"^ by 
modelling the global X-ray spectrum under the assumption that 
most of the emission is thermal. Spatially resolved spectral 
analysis showed that this assumption is only valid locally in 
the NW and in the SE. The global spectrum is in fact domi- 
nated by synchrotron emission. Therefore, even though their re- 
sult agrees with ours, we view this as coincidental. Numerical 
simulations have also been carried out to study the evolution 
of SNRs in constant density ISM. In order to be consistent 
with the observed angular size and expansion rate of SN 1006, 
Dwarkadas & Chevalier (1998) derived from their simulation a 
density between 0.05 cm"-' and 0. 1 cm"-'. Those results are com- 
patible with what we found. We note that our density is lower 
than values of most other studies, but most of them were fo- 
cused on the NW, where the ambient medium seems denser. The 
value found in the SE region is more likely to be representative 
of the overall ambient medium density, including in the bright 
synchrotron limbs. 

5.2. Limitations of our model 

Our method is based on estimating the density of the shocked 

ambient medium. To infer the density of the unshocked ambient 

medium, we have to assume a value for the compression ratio. 

In the efficient acceleration sites (NE and SW rims) the com- 

I ^ 1 

pressio n ratio can increase to 6 as calculated in Berezh ko et al.l 
(l2002h . However, in the SE and NW rims, the acceleration pro- 
cess is not as efficient, as we see very little X-ray or radio syn- 
chrotron emission. Thus a compression ratio of 4 seems reason- 
able. 

To derive the density from the EM, we assumed that the 
shocked ISM filled all the volume i.e. a filling factor e=l. As 
seen in Sect. 13.11 we need ejecta to reproduce the spectra and 
it must occupy some fraction of the volume. In the ejecta- 
dominated phase (Chevalier 1982), the fraction of the volume 
occupied by the shocked ISM in the global shocked region is 
about 0.8. SN 1006 is in a more advanced phase (towards Sedov) 
where the ejecta are more diluted in the shocked ISM. Moreover, 
closer to the shock, as in our regions, e is supposed to be closer 
to unity. So in our model, we can reasonably say that s = 0.8 is a 
lower limit to the filling factor. In that case, the density estimate 
for region ESE would increase from 0.053 cm"^ to 0.059 cm""*. 
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To derive the density from the norm parameter, we assumed a 
Sedov density profile to estimate the value of the integral J nidi. 
As SN 1006 is switching from the ejecta-dominated stage to the 
Sedov phase, the density profile that we assumed is steeper than 
it is in reality. We therefore overestimate the density correspond- 
ing to a given value of the EM. In the extreme case of a flat 
density profile, we obtain for the ESE region a density of 0.042 
cm"^ in comparison with 0.053 cm"^ for the Sedov profile. 

The di stance of 2.2 kpc us ed throughout this paper was 
derived by IWinkler et al.l ( l2003h from the measurement of the 
shock speed of the bright filament of 2890± 100 km/s and the 
appare nt expansion of 0".28 y r~'. However, lHeng & McCravl 
(12007'), using the same data as in lWinkler et all (l2003 l) but with a 
new model that incorporates additional physics, derived a lower 
shock speed of 2509±111 km/s. The corresponding distance to 
the remnant with this shock speed is 1 .9 kpc. For our study this 
will change the length of the line of sight in the EM method, but 
it will have a low impact on the density as «h ^ y/norm/dl. 

Concerning the contribution of the synchrotron emission, our 
three-component model was not able to disentangle the thermal 
emission of the shocked ISM from the synchrotron emission. We 
cannot exclude that part of the emission at high energy is due to 
synchrotron emission, so our density estimate is an upper limit 
in that respect. A lower limit to the density is discussed later in 
Sect. 123] 

5.3. Velocity of the blast wave 

With the density of 0.05 cm"^ inferred above, we can reestimate 
several parameters in order to account for the observed radius of 
9.6 pc (at a distance of 2.2 kpc). The age is t=1000 years, and 
as SN 1006 seems to be a type la supernova, the ejecta mass 
is Mej - 1.4 Mq. The impact of the distance uncertainty was 
low in the density estimate. However, parameters like energy or 
shock speed are really dependent on this distance. The following 
results are given for both d=2.2 kpc and 1 .9 kpc. As the remnant 
is just switching from the ejecta-dominate d stage to the Sedov 
stage, we decided to use the equation from iTruelove & McKe^ 
( 1 1 9991) for an homo geneous ISM with s=0 and a type la ejecta 
profile with n=7. 

The kinetic SN energy consistent with the previous pa- 
rameters is EsN=2. 2/1. 3x10^' ergs. For example, for a density 
higher than 0.2 cm the energy required would be higher than 
EsN=6.0/3.3 xlO^' ergs. Type la SN are expected to release an 
amount of energy of ~ 10^' ergs. Energies above 2x10^' ergs 
seem hard to reach. Thus a uniform gas density > 0.2 cm"^ sur- 
rounding all the remnant is excluded for energetic reasons. 

The global velocity of the blast wave can also be reestimated. 
As the density is lower in the SE, the shock speed there must 
be higher than the shock speed in the NW which is 2900 km/s. 
Indeed, the shock speed derived from the Truelove & McKee 
equations with an ambient medium density of 0.05 cm""* at 
t=1000 yeai-s is 4900/4600 km/s. If the density found in the SE 
region is representative of the rest of the remnant, this velocity 
is what we can expect for the whole remnant except in the NW. 
We note that there is also a difference between the global shock 
speed measured in radio and the velocity in the NW deduced 
from Ha observations. Indeed, even thou gh the error bars are 
big, the mean proper motion measured by lMoff"ett et al.l (Il993h 
in radio is 0".44±0".13 yr"', corresponding to 4500 km/s for a 
distance of 2.2 kpc. This is much higher than the velocity of 2890 
km/s obtained in the NW, and it is compatible with our estimate. 




I ^ ^ ^ I ^ ^ ^ I ^ ^ ^ I ^ ^ ^ I ^ ^ ^ 

0-CO 0.02 0.04 0.06 O.OB 0.10 

ambient medium density In part/cm-3 

Fig. 4. Theoretical and "model" ratio of pressures behind the 
transmitted shock in the NW cloud and behind the blast wave 
elsewhere. 



5.4. Lower limit on the density 

When a shock wave encounters a cloud, a transmitted shock is 
driven into the cloud. The pressure behind the transmitted shock 
is greater, at most, by a factor /3 than the pressure behind the 
blast wave not in interaction with the cloud. For a given density 
contrast ndoud/niSM, we can derive this /3 factor using the equa- 
tions from Sgro ( 1975). The pressure is proportional to nH-V^^^^^^- 
What we are interested in here are the density and velocity right 
at the shock. We know those values in the NW (transmitted 
shock) from the optical studies : »h - 0-25 - 0.4 cm~^ and 
KLck=2890km/s (Ravmon d et alj2007l:IWnkler et al.l200l re- 
spectively). For any ambient medium density, the velocity (in a 
region where the shock is not encountering any cloud) is known 
using the Truelove & McKee equations. Therefore we can de- 
rive a "model" ratio of pressures for any ambient medium den- 
sity. By comparing the theoretical ratio /3 with the "model" ratio, 
we can set a lower limit to the ambient medium density. The 
change of slope in the "model" ratio corresponds to the tran- 
sition between the ejecta-dominated phase and the Sedov phase. 
For n(;ioud=0.25 cm"-*, the "model" ratio becomes higher than the 
maximum theoretical value of /3(0.25/«ism) for nisM lower than 
0.03 cm"^ (see Fig. HJi. For ncioiid=0.4 cm"-*, the lower limit on 
the ambient medium density is 0.07 cm"-*. We conclude that the 
density of the ambient medium where the SNR evolves cannot 
be lower than 0.03 cm"^. For our preferred value of the density 
of 0.05 cm"-*, we have a "model" ratio of 1 .7 which is compatible 
with the theoretical value of j6(0.25/0.05)=2.1. 

5.5. Implication for CR acceleration 

The efficiency of the acceleration of the particles is strongly de- 
pendent on the shock speed. In SN 10 06, this velocity h a s been 
me asured only in the NW r egion, by iGhavamian et al.l (l2002h 
and iHeng & McCrayl ((2007). However, the velocity there does 
not seem to be representative of the rest of the remnant, and in 
particular, of the efficient acceleration sites in the NE and SW 
rims. If we use for those regions the velocity estimated in the SE 
in Sect. 15.31 this has an impact on the highest energy the accel- 
erator can reach. In fact, in a Bohm regime, at an age of 1000 
yr we have £,™x = ^^.Ij^ij^^f TeV (adapted from 
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iParizot et al.ll2006l) where is the downstream magnetic field. 
Thus with the higher velocity (4900 km/s) and Z?d=150 //G as 
in|Ksenofontov et al. (2005), the maximum energy that particles 
could reach in the bright limbs would be fimax - 1 PeV. 

The TeV emission of those accelerated particles depends on 
the density of the ambient medium. The total y-ray flux (IC+;7ro- 
decay) from the NE half of the remnant for different values of 
tin has been given in Ksenofontov et al. (2005). For an upstream 
magnetic field of 30 /uG, a density of 0.05 cm and a distance 
of d=2.2 kpc, this flux is a factor of 3 below the current HESS 
upper limit, and only a factor ~ 2 for d=1.9 kpc. This would, 
however, need to be revisited in the Ught of our revised estimate 
of the shock velocity. 

5.6. Electron temperature equilibration 

A question that remains open in the physics of strong shocks 
is the mechanism by which the electrons are heated just be- 
hind the shock. In the Rankine-Hugoniot relations for high 
Mach number shocks, the particles are heated in relation to 
their mass. The proton temperature is higher than the electron 
temperat ure by a factor 1836. However, it h as also been pro- 
posed by ICargill & PapadopoulosI (1 19881) and iGhavamian et alj 
( l2007h that collective plasma interaction can give rise to prompt 
(limited) electron heating at the shock. Such mechanisms pro- 
duce a higher electron temperature at the shock than in the 
Rankine-Hugoniot case. Our model requires an electron temper- 
ature kT^e ~ 2 keV to reproduce the spectrum in the region ESE. 
The Rankine-Hugoniot relations (without particle acceleration) 
imply a proton temperature just at the shock of ~ 40 keV for a 
shock speed of 4900 km/s, if the degree of temperature equi- 
libration is low. Using the Coulomb equilibration mechanism 
with an initial electron-to-proton temperature ratio of 1/1836 
(the electron-to-proton mass ratio), the predicted electronic tem- 
perature is lower than what is observed. Indeed, with an ioniza- 
tion timescale of 5 x 10^ cm"^ s, kT^ rises to 0.8 keV only. An ini- 
tial ratio of at least re/ri=0.04 (corresponding to kre=1.7 keV 
just at the shock) is needed for the temperature to reach 2 keV. 
According to Ghavamian et al. (2007), the electronic tempera- 
ture immediately behind the shock is constant (~ 0.3 keV) for 
any shock speed above 400 km/s. In comparison, the electronic 
temperature just at the shock required in our model (kT^e^l.? 
keV) is much higher. 

However, we note that this high temperature depends very 
much on the amount of synchrotron emission in the SE, as 
discussed in Sect. 13.51 To check if the three-component (syn- 
chrotron, shocked ejecta and shocked ISM) model is compat- 
ible with a low electronic temperature just at the shock, we 
fixed several parameters of the shocked ISM component for the 
lowest possible value of the ambient medium density i.e. 0.03 
cm"^. The norm parameter and the upper limit of the ioniza- 
tion timescale can be fixed from Sedov models and from the 
geometry of the emitting region. They are fixed respectively to 
0.8x10"^ cm"^ and 3 x 10^ cm""* s. From the upper limit of 
T, we can derive an upper limit to the temperature but not the 
average (for the region) electron temperature that is needed for 
the PSHOCK model. In order to avoid this difficulty, we use a 
plane-parallel shock model with separate ion and electron tem- 
peratures (NPSHOCK) to reproduce the shocked ISM. In this 
model the two temperatures can be set precisely. The mean post- 
shock temperature is derived from the shock speed (22 keV at 
4900 km/s), and the electron temperature immediately behind 
the shock front is set to 0.3 keV. The global fit is good and the 
C-stat is higher only by AC-stat=5 in comparison to the three- 



component model (where the parameters for the shocked ISM 
component were left free). Thus the data are compatible with a 
low electron temperature just behind the shock if a model in- 
cluding synchrotron emission is assumed. 

6. Summary 

The study of the post-shock thermal emission in the SE and in 
the NW of SN 1006 with XMM-Newton observations leads to the 
following conclusions : 

1 . Most previous studies were focused on the NW region where 
the remnant seems to be encountering a denser ISM. By 
modelling the post-shock thermal X-ray emission, we were 
able to estimate the density in the SE rim. The low value of 
~ 0.05 cm""* we found seems representative of the rest of the 
remnant. 

2. The shock speed in the SE region, and maybe in the rest of 
the remnant, is ~ 4900 km/s. The velocity of 2890 km/s was 
deduced from UV observations in the NW where the density 
is higher, causing the expansion to slow down. 

3. The highest energy particles can reach scales as V^^^^^. The 
new estimate of the velocity could increase E^ax up to ener- 
gies of about 1 PeV. 

4. If we assume that the ambient medium density found in the 
SE region extends to the brigh t limbs where the accelera tion 
is efficient, then according to iKsenofontov et al.] (l2005h the 
TeV gamma-ray flux would be a factor of ~ 3 below the cur- 
rent HESS upper limit. 

Acknowledgements. We are grateful to the referee for his valuable suggestions 
which helped us to improve the clarity of the paper 



References 

Ahai-onian, R, Akhperjanian, A., Banio, J., et al. 2001, A&A, 370, 112 
Aharonian, R, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2005, A&A, 437, 
L7 

Aharonian, R, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 449, 
223 

Berezhko, E. G., Ksenofontov, L. T., & Volk, H. J. 2002, A&A, 395, 943 
Cargill, R J. & Papadopoulos, K. 1988, ApJ, 329, L29 
Carter, J. A. & Read, A. M. 2007, A&A, 464, 1 155 
Chevalier R. A. 1982, ApJ, 258, 790 

Dubner, G. M., Giacani, E. B., Goss, W. M., Green, A. J., & Nyman, L. A. 2002, 

A&A, 387, 1047 
Dwarkadas, V. V. & Chevalier, R. A. 1998, ApJ, 497, 807 
Dyer, K. K., Reynolds, S. R, Borkowski, K. J., Allen, G. E., & Petre, R. 2001, 

ApJ, 551,439 

Ghavamian, R, Laming, J. M., & Rakowski, C. E. 2007, ApJ, 654, L69 
Ghavamian, R, Winkler R R, Raymond, J. C, & Long, K. S. 2002, ApJ, 572, 
888 

Hamilton, A. J. S., Sarazin, C. L., & Szymkowiak, A. E. 1986, ApJ, 300, 698 
Heng, K. & McCray, R. 2007, ApJ, 654, 923 

Iwamoto, K., Brachwitz, R, Nomoto, K., et al. 1999, ApJS, 125, 439 
Korreck, K. E., Raymond, J. C, Zurbuchen, T. H., & Ghavamian, R 2004, ApJ, 
615, 280 

Koyama, K., Petre, R., Gotthelf, E. V., et al. 1995, Nature, 378, 255 
Ksenofontov, L. T., Berezhko, E. G., & Volk, H. J. 2005, A&A, 443, 973 
Laming, J. M., Raymond, J. C, McLaughlin, B. M., & Blair, W. R 1996, ApJ, 

472, 267 
Lodders, K. 2003, ApJ, 591, 1220 

Long, K. S., Reynolds, S. P, Raymond, J. C, et al. 2003, ApJ, 586, 1 162 
Molfett, D. A., Goss, W. M., & Reynolds, S. P 1993, AJ, 106, 1566 
Parizot, E., Marcowith, A., Ballet, J., & Gallant, Y. A. 2006, A&A, 453, 387 
Raymond, J. C, Korreck, K. E., Sedlacek, Q. C, et al. 2007, ApJ, 659, 1257 
Reynolds, S. P & Gilmore, D. M. 1986, AJ, 92, 1138 
Rothenflug, R., Ballet, J., Dubner, G., et al. 2004, A&A, 425, 121 
Sgro, A. 1975, ApJ, 197, 621 

Tmelove, J. K. & McKee, C. R 1999, ApJS, 120, 299 
Wang, C. Y. & Chevalier, R. A. 2001, ApJ, 549, 1119 



F. Acero, J. Ballet, A. Decourchelle: The gas density around SN 1006 



9 



Winkler, P. R, Gupta, G., & Long, K. S. 2003, ApJ, 585, 324 
Winkler, P. F. & Long, K. S. 1997, ApJ, 491, 829 



